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Abstract 

We present a systematic study of the properties of pure hadronic and hybrid 
compact stars. The nuclear equation of state (EoS) for /3-equilibrated neutron star 
matter was obtained using density dependent effective nucleon-nucleon interaction 
which satisfies the constraints from the observed flow data from heavy-ion colli- 
sions. The energy density of quark matter is lower than that of this nuclear EoS 
at higher densities implying the possibility of transition to quark matter inside the 
core. We solve the Einstein's equations for rotating stars using pure nuclear matter 
and quark core. The /3- equilibrated neutron star matter with a thin crust is able 
to describe highly massive compact stars but find that the nuclear to quark mat- 
ter deconfinement transition inside neutron stars causes reduction in their masses. 
Recent observations of the binary millisecond pulsar J1614-2230 by P. B. Demorest 
et al. [1] suggest that the masses lie within 1.97=b0.04 Mq where Mq is the solar 
mass. In conformity with recent observations, pure nucleonic EoS determines that 
the maximum mass of NS rotating with frequency below r-mode instability is ~1.95 
Mq with radius ~10 kilometers. Although compact stars with quark cores rotating 
with Kepler's frequency have masses up to ~2 Mq, but if the maximum frequency 
is limited by the r-mode instability, the maximum mass ~1.7 Mq turns out to be 
lower than the observed mass of 1.97ib0.04 Mq, by far the highest yet measured 
with such certainty, implying exclusion of quark cores for such massive pulsars. 

PACS numbers: 26.60.-c, 21.65. Cd, 21.65.Ef, 26.60.Kp, 12.38.-t, 12.39.-x, 21.65. Qr 
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1 Introduction 



The mass-radius relationship of neutron stars (NSs) is of prime importance to 
understand the high-density low-temperature regime of the hadronic equation 
of state (EoS). Depending on this relationship, certain models for the hadronic 
EoS can either be confirmed or ruled out. Several attempts have been made 
on measuring the radii and masses of NSs to constrain the uncertainties in 
the high density behavior of the EoS. The observations on double NSs [2], 
ghtches in radio pulsars [3] , thermal emission [4] from accreting NSs and from 
millisecond X-ray pulsars lead to constraints on mass-radius relationship of 
NSs. Recently the pressure of neutron star matter at supranuclear density 
is measured by Ozel et al. [5] directly from observations using advanced as- 
trophysical techniques and NS atmosphere modeling. The pressure extracted 
from NS mass-radius data crucially constrains the extension of the EoS to high 
density low temperature regime for stellar matter ruling out those who fail to 
comply with the recent mass-radius data. The quark matter can support a 
star as massive as ~2 Mq only if the quarks are strongly interacting and are 
therefore not 'free' quarks. To overcome this situation, Dexheimer et al. [6,7] 
have recently employed a hadronic SU(3) sigma-omega model including Delta- 
resonances and hyperons to describe the properties of neutron star matter by 
softer equation of state. Delta-resonances have a repulsive vector potential 
which works to counteract gravity in a compact star. They successfully repro- 
duce both the measured mass-radius relationship and the extrapolated EoS by 
slightly lowering the coupling strength of the Delta resonances to the vector 
mesons. 



In our previous works, the density dependent M3Y effective interaction (DDM3Y) 
which provides a unified description of the elastic and inelastic scattering [8,9], 
cluster [10], a [10,11,12] and proton radioactivities [13], the symmetric and 
asymmetric nuclear matter [13,14,15], was employed to obtain nucleonic EoS 
of the /3-equihbrated NS matter [16,17]. At high densities, the energy density 
of this ^-equilibrated charge neutral NS matter is higher than that of quark 
matter signalling the onset of deconfinement transition to quark matter inside 
the star. In the present work, we solve the Einstein's equations for rotating 
stars using pure nuclear matter without and with quark matter core. A sys- 
tematic study of the static as well as rotating compact stars with quark matter 
inside is presented in view of the recent observations of the massive compact 
stars. We shall find later that the present EoS unlike other EoS [18,19,20,21] 
can explain successfully the recently observed mass-radius data. The effect of 
the presence of the quark core on the determination of maximum mass of NS 
will be investigated for both static and rotating stars. 
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2 Present observational status 



With the energies and interaction rates foreseen at FAIR, the compressed 
baryonic matter (CBM) will create highest baryon densities in nucleus-nucleus 
collisions to explore the properties of superdense baryonic matter and the in- 
medium modifications of hadrons. The compact stars provide natural testing 
laboratory for highly compressed matter. The stiffness of the high-density 
matter controls the maximum mass of compact stars. The analyses of mass- 
radius data on NSs by Ozel et al. [5] favor smaller masses lying within 1.6- 
1.9 M0 with radii 8-10 kilometers. Recent mass measurement of the binary 
millisecond pulsar J1614-2230 by P.B. Demorest et al. [1] rules out the EoS 
which fail to predict the masses within 1.97±0.04 Mq. Most of the currently 
proposed EoS [18,19,20,21] involving exotic matter, such as kaon condensates 
or hyperons failed to produce such a massive star. The measured mass of 
PSR J1748-2021B, a millisecond pulsar in the Globular Cluster NGC 6440, is 
claimed to be as high as 2.74lo;|];M0 (2cr) [22]. Moreover, a pulsar rotating 
faster (e.g., PSR J17482446ad) than the limit set by the r-mode instability has 
already been observed [23] . Further observations and a better r-mode modeling 
may shed more light on this issue. 



3 Construction of compact star models 

If rapidly rotating compact stars were nonaxisymmetric, they would emit 
gravitational waves in a very short time scale and settle down to axisymmet- 
ric configurations. Therefore, we need to solve for rotating and axisymmetric 
configurations in the framework of general relativity. For the matter and the 
spacetime the following assumptions are made. The matter distribution and 
the spacetime are axisymmetric, the matter and the spacetime are in a sta- 
tionary state, the matter has no meridional motions, the only motion of the 
matter is a circular one that is represented by the angular velocity, the angular 
velocity is constant as seen by a distant observer at rest and the matter can be 
described as a perfect fluid. The energy-momentum tensor of a perfect fluid 
T^'^ is given by 

T^"" = {e + P)u^'vr -g^T (1) 

where £, P, and g'^^ are the energy density, pressure, four velocity and the 
metric tensor, respectively. To study the rotating stars the following metric is 
used 

^e^i-p)r'^ sin2 Qi^dcf) - udtf (2) 
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where the gravitational potentials 7, p, a and u are functions of polar coor- 
dinates r and 6 only. The Einstein's field equations for the three potentials 7, 
p and a have been solved using the Green 's-f unction technique [24] and the 
fourth potential co has been determined from other potentials. All the phys- 
ical quantities may then be determined from these potentials. Obviously, at 
the zero frequency limit corresponding to the static solutions of the Einstein's 
field equations for spheres of fluid, the present formalism yields the results for 
the solution of the Tolman-Oppenheimer-Volkoff (TOV) equation [25]. We use 
the 'rns' code [26] for calculating the compact star properties which requires 
EoS in the form of energy density versus pressure along with corresponding 
enthalpy and baryon number density and since we are using various EoS for 
different regions, these are smoothly joined. 



4 The crustal equations of state 

The different regions of a compact star are governed by different EoS. These 
can be broadly divided into two regions: a crust that accounts for about 5% 
of mass and about 10% of the radius of a star and the core is responsible for 
the rest of the mass and radius of a star. The outer layers are a solid crust 
~ 1 km thick, consisting, except in the outer few meters, of a lattice of bare 
nuclei immersed in a degenerate electron gas. As one goes deeper into the 
crust, the nuclear species become, because of the rising electron Fermi energy, 
progressively more neutron rich, beginning (ideally) as ^^Fe through ^-"^^Kr at 
mass density f«4.3xl0^^ g cm~^. At this density, the 'neutron drip' point, the 
nuclei have become so neutron rich that with increasing density the continuum 
neutron states begin to be filled, and the lattice of neutron- rich nuclei becomes 
permeated by a sea of neutrons. 

The EoS that cover the crustal region of a compact star are Feynman-Metropolis- 
Teller (FMT) [27], Baym-Pethick-Sutherland (BPS) [28] and Baym-Bethe- 
Pethick (BBP) [29]. The most energetically favorable nucleus at low densities 
is ^^Fe, the endpoint of thermonuclear burning. The FMT is based on Fermi- 
Thomas model to derive the EoS of matter at high pressures and covers the 
outermost crust which is essentially made up of iron and a fraction of the 
electrons bound to the nuclei. The major difficulty in deriving the equation of 
state is the calculation of the electronic energy. At subnuclear densities, from 
about 10"^ g cm~'^ up to the neutron drip density 4.3x10^^ g cm~^ the EoS 
of BPS is apphcable which includes the effects of the lattice Coulomb energy 
on the equilibrium nuclide. The domain from neutron drip density to about 
nuclear density 2.5x10^^ g cm~^, is composed of nuclei, electrons and free 
neutrons where EoS of BBP is applicable which is based on a compressible 
liquid drop model of nuclei with conditions that nuclei must be stable against 
^-decay and free neutron gas must be in equilibrium with neutrons in nuclei. 
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5 The ^-equilibrated nucleEir matter and quark matter EoS for 
compact star core 



The nuclear matter EoS is calculated [13] using the isoscalar and the isovector 
components of M3Y interaction along with density dependence. The density 
dependence of the effective interaction, DDM3Y, is completely determined 
from nuclear matter calculations. The equilibrium density of the nuclear mat- 
ter is determined by minimizing the energy per nucleon. The energy variation 
of the zero range potential is treated accurately by allowing it to vary freely 
with the kinetic energy part e'^™ of the energy per nucleon e over the entire 
range of e. In a Fermi gas model of interacting neutrons and protons, the 
energy per nucleon for isospin asymmetric nuclear matter [13] is given by 

<P.X) = [^]F[X) + [P^)[l-Ppl) (3) 



where kp={l.^Ti'^ p)^ which equals Fermi momentum in case of symmetric nu- 
clear matter (SNM), the isospin asymmetry X = with p = Pn + Pp, Pn 
and Pp being the neutron and proton number densities respectively, and the 

kinetic energy per nucleon e'=^"=[^^]F(X) with F(X) = [^i±^^^^^±^i^=^]. 
Jv—Jvoo + ^^Jvoi, Jvoo and J^oi represent the volume integrals of the isoscalar 
and the isovector parts of the M3Y interaction. The isoscalar t^^^ and the 
isovector t^^^ components of M3Y interaction potential are, respectively, 
given by t^^y^^^^fcm^ ^ 799922^ _ 21342st|5£) + j^^^i _ ae^-)5(s) 
and t^f^y^g^g/cm) ^ _488Q22Pl=M ^ n^Q^^^st^ '+ Joi{l - ae'''')S{s) where 
Joo = -276 MeVfm3, Jqi = 228 MeVfm^, a = 0.005 MeV-^ The DDM3Y 
effective NN interaction is given by voi{s,p,e) — t^^^ {s,€)g{p) where the 
density dependence g{p) = C(l — /Jpa) and the constants C and /3 of den- 
sity dependence have been obtained from the saturation condition |^ = at 
p = Po and e = eo where po and eo are the saturation density and the satu- 
ration energy per nucleon respectively [13]. This EoS evaluated at the isospin 
asymmetry X determined from the /5-equilibrium proton fraction [— y], 

obtained by solving hc{37r'^pxpy^^ — 4£'sj,^(p)(l — 2x13), provides EoS for the 
/3-equilibrated NS matter where nuclear symmetry energy Egymip) is given by 

i?.,m(p) = e(p,l)-e(p,0). (4) 



The areas enclosed by the continuous and dashed hnes in Fig.-l correspond to 
the pressure regions for neutron matter consistent with the experimental flow 
data after inclusion of the pressures from asymmetry terms with weak (soft 
NM) and strong (stiff NM) density dependences, respectively [30]. Although, 
the parameters of the density dependence of DDM3Y interaction are tuned to 
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1 2 3 4 5 6 7 

P/Po 



Fig. 1. The pressures P of the /3-equihbrated charge- neutral neutron star matter 
[DDM3Y(/3eq)] and pure neutron matter for the present work [DDM3Y] and that of 
Akmal et al. [31] are plotted as functions of density p/ Pq. The areas enclosed by the 
continuous and dashed lines correspond to the pressure regions for neutron matter 
consistent with the experimental flow data after inclusion of the pressures from 
asymmetry terms with weak (soft NM) and strong (stiff NM) density dependences, 
respectively [30]. 



reproduce the saturation energy per nucleon eo and the saturation density po 
of symmetric nuclear matter which are obtained from finite nuclei, the agree- 
ment of this EoS with the experimental flow data [30], where the high density 
behaviour looks phenomenologically confirmed, justifies its extrapolation to 
high density [13]. 



For cold and dense quark (QCD) matter, the perturbative EoS [32] with two 
massless and one massive quark flavors and a running coupling constant, is 
used. The constant B is treated as a free parameter, which allows to take into 
account non-perturbative effects not captured by the weak coupling expansion. 
In fact, using the free quark number density, one recovers the expression for 
the pressure in the original MIT bag model [33] , with B taking the role of the 
bag constant. Due to physics criteria (e.g. requiring the energy density to be 
positive), the possible values for B are, however, typically rather restricted, 
allowing to make quantitative statements that are not possible in the original 
MIT bag model. 
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Fig. 2. Variation of mass with central density for static and rotating neutron stars 
with pure nuclear matter inside. 
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Fig. 3. Mass-radius relationship for static and rotating neutron stars with pure 
nuclear matter inside. 



6 Calculations and Results 



The rotating compact star calculations are performed using the crustal EoS, 
FMT + BPS + BBP upto number density of 0.0458 fm~^ and /3-equilibrated 
NS matter beyond. It is worthwhile to mention here that a star may not 
rotate as fast as Keplerian frequency due to r-mode instability. There have 
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Fig. 4. The EoS of the /3-equihbrated charge neutral neutron star matter and the 
quark matter EoS. 

been suggestions that the r-mode instability may Umit the time period to 1.5 
ms [34]. However, a pulsar rotating faster (e.g., PSR J17482446ad) than this 
limit has already been observed [23] . The variation of mass with central density 
for static and rotating neutron stars at Keplerian limit and also maximum 
frequencies limited by the r-mode instability with pure nuclear matter inside is 
shown in Fig.-2. In Fig.-3, the mass-radius relationship for static and rotating 
neutron stars at Keplerian limit and also at maximum frequencies limited by 
the r-mode instability with pure nuclear matter inside is shown. Fig.-3 depicts 
that NSs with pure nuclear matter inside, the maximum mass for the static 
case is 1.92 with radius ~9.7 km and for the star rotating with Kepler's 
frequency it is 2.27 Mq with equatorial radius ~13.1 km [16]. However, for 
stars rotating with maximum frequency limited by the r-mode instability, the 
maximum mass turns out to be 1.95 (1.94) corresponding to rotational 
period of 1.5 (2.0) ms with radius about 9.9 (9.8) kilometers [17]. 

The energy density of the quark matter is lower than that of the present 
EoS for the /^-equilibrated charge neutral NS matter at densities higher than 
0.405 fm~^ for bag constant i?4=110 MeV [32] implying presence of quark 
core. The energy densities of the present EoS for the /^-equilibrated charge 
neutral NS matter and the quark matter EoS for bag constant i?4=110 MeV 
are shown in Fig.-4 as functions of baryonic densities. For lower values of 
bag constant such as i?4=89 MeV, energy density for our EoS is lower and 
makes a cross over with the quark matter EoS at very high density ~1.2 
fm~^ causing too little quark core and therefore we choose i? 2=110 MeV 
for representative calculations. The common tangent is drawn for the energy 
density versus density plots where pressure is the negative intercept of the 
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Fig. 5. Variation of mass with central density for static and rotating neutron stars 
with nuclear and quark matter inside. 

tangent to energy density versus density plot. However, as obvious from Fig.- 
4, the phase co-existence region is neghgibly small which is represented by part 
of the common tangent between the points of contact on the two plots [35] 
implying constant pressure throughout the phase transition. The variation of 
mass with central density for static and rotating neutron stars at Keplerian 
limit and also maximum frequencies limited by the r-mode instability with 
nuclear and quark matter inside is shown in Fig.-5. In Fig.-6, the mass-radius 
relationship for static and rotating neutron stars at Keplerian limit and also at 
maximum frequencies limited by the r-mode instability with nuclear and quark 
matter inside is shown. Fig.-6 depicts that when quark core is considered, the 
maximum mass for the static case is 1.68 Mq with radius ~10.4 km and for 
the star rotating with Kepler's frequency it is 2.02 M© with equatorial radius 
~14.3 km. In a similar study, it was concluded that compact stars with a 
quark matter core and an hadronic outer layer, can be as massive as 2.0 M© 
but stay below the pure quark stars and pure neutron stars [36]. However, 
they have used two different relativistic mean-field parameter sets TMl and 
NL3 [36] to explore the influence of the hadronic part of the EoS whose high 
density behaviour do not satisfy the criteria extracted from the experimental 
flow data [30]. For our case stars rotating with maximum frequency limited 
by the r-mode instability, the maximum mass turns out to be 1.72 (1.71) M© 
corresponding to rotational period of 1.5 (2.0) ms with radius about 10.7 (10.6) 
kilometers. 

It is worthwhile to mention here that the consideration of rotation at the 
Kepler frequency is not of much relevance for a comparison with observational 
data since even the shortest observed rotation period is still far from the 
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Fig. 6. Mass-radius relationship for static and rotating neutron stars with nuclear 
and quark matter inside. 

Keplerian limit and does not result in a sufficient increase in the maximum 
mass. The recent observation of PSR J1614-2230 with a mass of 1.97±0.04 
Mq rotating with a period of 3.1 ms shows that either an EoS or transition to 
quark matter should be excluded that does not allow for star configurations 
with maximum mass reaching such high values. Since the hadronic EoS used 
in the present work is not soft (incompressibility 274.7±7.4 MeV) [13] at high 
densities and neutron star masses reach as high as ~2 Mq for such rotational 
frequencies, implies that theoretically robust quark matter EoS [32] used in 
this work is not stiff enough at high densities which was obtained from ffist 
principle calculations based on perturbation theory taking terms up to 0{a1) 
with quark chemical potentials and strange quark mass non-zero. The effects of 
quark pairing were incorporated by adding to the pressure a term accounting 
for the condensation energy of the cooper pairs in the color-flavor-locked phase. 
However, a more detailed and realistic inclusion of the effect of pairing may 
make quark EoS stiffer. 



7 Summary and Conclusion 



In summary, the energy density of the present EoS for /3-equilibrated charge 
neutral NS matter using DDM3Y effective NN interaction turns out to be 
higher than that of quark matter at densities above 0.405 fm~^ implying pos- 
sibility of quark core. We have applied our nucleonic EoS with a thin crust 
to solve the Einstein's field equations to determine the mass-radius relation- 
ship of neutron stars with and without quark cores. We have obtained the 
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masses of neutron (hybrid) stars rotating with Keplerian frequencies, around 
2.27 (2.02) Mq with equatorial radii around 13 (14) kilometres. The result 
for NS without quark core is in excellent agreement with recent astrophysi- 
cal observations. The neutron star matter can further undergo deconfinement 
transition to quark matter, thereby reducing compact star masses consider- 
ably. Although such hybrid compact stars rotating with Kepler's frequency 
have masses up to ~2 M0, but rotating with maximum frequency limited by 
the r-mode instability (or with 3.1 ms as observed for pulsar J1614-2230), the 
maximum mass ~1.7 turns out to be lower than the observed mass of 
1.97±0.04 Mq and thus rules out quark cores for such massive pulsars but 
not for pulsars with masses ~1.7 Mq or less. Obviously, in order not to con- 
flict with the mass measurement [1], either there must be some mechanism 
to prevent nuclear matter to deconfine into quark matter or the quark EoS 
should be made stiffer by several possible realistic improvements [32]. The 
nucleon-nucleon effective interaction used in the present work, which is found 
to provide a unified description of elastic and inelastic scattering, various ra- 
dioactivities and nuclear matter properties, also provides an excellent descrip- 
tion of the /3-equilibrated NS matter which is stiff enough at high densities to 
reconcile with the recent observations of the massive compact stars ~2 Mq 
while the corresponding symmetry energy is supersoft [14,15] as preferred by 
the FOPI/GSI experimental data. 
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